Precision Measurement of Neutrino Oscillation Parameters with KamLAND 
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The KamLAND experiment has determined a precise value for the neutrino oscillation parameter Am|i and 
stringent constraints on 9\2- The exposure to nuclear reactor anti-neutrinos is increased almost fourfold over pre- 
vious results to 2.44x 10'^^ proton-yr due to longer livetime and an enlarged fiducial volume. An undistorted re- 
actor Ue energy spectrum is now rejected at >5ct. Analysis of the reactor spectrum above the inverse beta decay 
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energy treshold, and including geo-neutrinos, gives a best-fit at Am2i = 7.58 
and tan2 6li2=0.56i;5 o? (stat)lQ Qg(syst). Local Ax^-minima at higher and lower Am2i are disfavored at 
>4a. Combining with solar neutrino data, we obtain Amii =7.59lg;2i x 10" 



'eV^ and tan^6li2 = 0.47: 
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PACS numbers: 14.60.Pq, 26.65.+t, 28.50.Hw, 91.35.-x 

Experiments studying atmospheric, solar, reactor and ac- 
celerator neutrinos provide compelling evidence for neutrino 
mass and oscillation. The Kamioka Liquid scintillator Anti- 
Neutrino Detector (KamLAND) investigates neutrino oscil- 
lation parameters by observing electron anti-neutrinos (Fg) 
emitted from distant nuclear reactors. Previously, KamLAND 
announced the first evidence of Vf. disappearance [1], fol- 
lowed by direct evidence for neutrino oscillation by observing 
distortion of the reactor 77^ energy spectrum [2]. More re- 
cently, KamLAND showed the first indication of geologically 
produced anti-neutrinos (geo-neutrinos) from radioactive de- 
cay in the Earth [3], possibly a unique tool for geology. 

This Letter presents a precise measurement of Amji and 
new constraints on 9i2 based on data collected from March 
9, 2002 to May 12, 2007, including data used earlier [1,2]. 
We have enlarged the fiducial volume radius from 5.5 m to 



6 m and collected significantly more data; the total exposure 
is 2.44x lO^^^ proton-yr (2881 ton-yr). We have expanded the 
analysis to the full reactor energy spectrum and reduced 
the systematic uncertainties in the number of target protons 
and the background. We now observe almost two complete 
oscillation cycles in the spectrum and extract more precise 
values of the oscillation parameters. 

KamLAND is at the site of the former Kamiokande experi- 
ment at a depth of ^^2700 m water equivalent. The heart of the 
detector is 1 kton of highly purified liquid scintillator (LS) en- 
closed in an EVOH/nylon balloon suspended in purified min- 
eral oil. The LS consists of 80% dodecane, 20% pseudoc- 
umene and 1.36 ± 0.03 g/1 of PPO [4]. The anti-neutrino de- 
tector is inside an 18-m-diameter stainless steel sphere. An 
array of 1879 50-cm-diameter photomultiplier tubes (PMTs) 
is mounted on the inner surface of the sphere. 554 of these are 
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TABLE I: Estimated systematic uncertainties relevant for the neu- 
trino oscillation parameters Am|i and 612. 





Detector-related (%) 




Reactor-related (%) 






Energy scale 


1.9 


77e-spectra [7] 


0.6 




Fiducial volume 


1.8 


I7e-spectra 


2.4 


Event rate 


Energy threshold 


1.5 


Reactor power 


2.1 




Efficiency 


0.6 


Fuel composition 


1.0 




Cross section 


0.2 


Long-lived nuclei 


0.3 



reused from the Kamiokande experiment, while the remain- 
ing 1325 are a faster version masked to 17 inches. A 3.2-kton 
cylindrical water-Cherenkov outer detector (OD), surrounding 
the containment sphere, provides shielding and operates as an 
active cosmic-ray veto detector. 

Electron anti-neutrinos are detected via inverse /3-decay, 
Ve+p e'^ + 71, with a 1 .8 MeV threshold. The prompt scin- 
tillation light from the e"*" gives a measure of the Ve energy, 
~ E'p + En + 0.8 MeV, where E^, is the prompt event 
energy including the positron kinetic and annihilation energy, 
and En is the average neutron recoil energy, O(lOkeV). The 
mean neutron capture time is 207.5 ± 2.8 ^s. More than 99% 
capture on free protons, producing a 2.2 MeV 7 ray. 

KamLAND is surrounded by 55 Japanese nuclear power 
reactor units, each an isotropic source. The reactor oper- 
ation records, including thermal power generation, fuel bur- 
nup, and exchange and enrichment logs, are provided by 
a consortium of Japanese electric power companies. This 
information, combined with publicly available world re- 
actor data, is used to calculate the instantaneous fission 
rates using a reactor model [5]. Only four isotopes con- 
tribute significantly to the Ve spectra; the ratios of the fis- 
sion yields averaged over the entire data taking period are: 
235u : 238u : 239pu ; 24ipy ^ o.570: 0.078: 0.295: 0.057. The 
emitted Ve energy spectrum is calculated using the Ve spectra 
inferred from Ref. [6], while the spectral uncertainty is eval- 
uated from Ref. [7]. We also include contributions from the 
long-lived fission daughters ^"Sr, ^"®Ru, and ^''^Ce [8]. 

We recently commissioned an "off-axis" calibration sys- 
tem capable of positioning radioactive sources away from the 
central vertical axis of the detector. The measurements in- 
dicate that the vertex reconstruction systematic deviations are 
radius- and zenith-angle-dependent, but smaller than 3 cm and 
independent of azimuthal angle. The fiducial volume (FV) is 
known to 1.6% uncertainty up to 5.5 m using the off-axis cal- 
ibration system. The position distribution of the /3-decays of 
muon-induced ^^B/^^N confirms this with 4.0% uncertainty 
by comparing the number of events inside 5.5 m to the num- 
ber produced in the full LS volume. The ^^B/^^N event ratio 
is used to establish the uncertainty between 5.5 m and 6m, 
resulting in a combined 6-m-radius FV uncertainty of 1.8%. 

Off-axis calibration measurements and numerous central- 
axis deployments of ^OCo, '^^Ge, ^"^^Hg, ^'^Zn, ^'^'^Am^Be, 
^^^Cs and ^^"Po^^C radioactive sources established the event 



reconstruction performance. The vertex reconstruction res- 
olution is ~12cm/^/E(MeV) and the energy resolution is 
6.5%/^E{MeY). The scintillator response is corrected for 
the non-linear effects from quenching and Cherenkov light 
production. The systematic variation of the energy recon- 
struction over the data-set give an absolute energy-scale un- 
certainty of 1.4%; the distortion of the energy scale results 
in a 1.9% uncertainty on Am^i, while the uncertainty at the 
analysis threshold gives a 1.5% uncertainty on the event rate. 
Table I summarizes the systematic uncertainties. The total 
uncertainty on Am|i is 2.0%, while the uncertainty on the 
expected event rate, which primarily affects 9i2, is 4.1%. 

For the analysis we require 0.9 MeV< Ep < 8.5 MeV. The 
delayed energy, E^, must satisfy 1.8 MeV < E^ < 2.6 MeV 
or 4.0 MeV < Ed < 5.8 MeV, corresponding to the neutron- 
capture 7 energies for p and ^^C, respectively. The time 
difference (AT) and distance (Ai?) between the prompt 
event and delayed neutron capture are selected to be 0.5 /xs< 
AT < 1000 /is and Ai? < 2 m. The prompt and delayed radiiil 
distance from the detector center {Rp, Ra) must be <6 m. 

Accidental coincidences increase near the balloon surface 
(i? = 6.5m), reducing the signal-to-background ratio. We 
use constraints on event characteristics to suppress acciden- 
tal backgrounds while maintaining high efficiency. We con- 
struct a probabiUty density function (PDF) for accidental co- 
incidence events, facc{Ep, E^, AR, AT, Rp, R^), by pairing 
events in a lO-ms-to-20-s delayed-coincidence window. A 
PDF for the signal, fj^^ {Ep, E^, AR, AT, Rp, Ri), is con- 
structed from a Monte Carlo simulation of the prompt and 
delayed events using the measured neutron capture time and 
detector response. For the Ep distribution in /i7^ , we choose 
an oscillation-free reactor spectrum including a contribution 
from geo-neutrinos estimated from Ref. [9]. A discriminator 
value, L = -jrr^- — . is calculated for each candidate pair 
that passes the earlier cuts. We establish a selection value 
j^cut i^jj^g 0.1 MeV, where L^"* is the value of L at 
which the figure-of-merit, ^g.i^g. is maximal. Si is the num- 
ber of Monte Carlo signal events in the ith energy bin with 
L > L?"*. Bi is calculated similarly using the accidental 
coincidence event pairs. The choice of the Ep distribution in 
/i7_. affects only the discrimination power of the procedure; 
substituting the oscillation-free reactor spectrum by an oscil- 
lated spectrum with the parameters from Ref. [2J changes our 
oscillation parameter results by less than Q.la. The selection 
efficiency e(£'p) is estimated from the fraction of selected co- 
incidence events relative to the total generated in R < 6m in 
the simulation, see Fig. l(top). 

The dominant background is caused by ^"^C(Q;,n)^^0 reac- 
tions from a-decay of ^^°Po, a daughter of -^^^Rn introduced 
into the LS during construction. We estimate that there are 
(5.56 ± 0.22) X 10'' 2iOpo a-decays. The ^^C(a,ny^O re- 
action results in neutrons with energies up to 7.3 MeV, but 
most of the scintillation energy spectrum is quenched be- 
low 2.7 MeV. In addition, i2c(n,n')^^C*, and the P' and 2"'* 
excited states of ^®0 produce signals in coincidence with 
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TABLE II: Estimated backgrounds after selection efficiencies. 



Selection efficiency 



Background 



Contribution 



Accidentals 
^Li/*He 

Fast neutron & Atmospheric v 
"C(a,n)i«05, 
"C(Q,n)i''09, 



80.5 ±0.1 
13.6±1.0 
<9.0 

np^np 157.2 ±17.3 

^^C(n,n')^2C* (4.4MeV7) 6.1 ±0.7 



^C(a,n)^*'0 P' exc. state (6.05 MeV e+e") 15.2 ± 3.5 
'C(Q,n)i'^0 2"'' exc. state (6. 13 MeV 7) 3.5 ±0.2 



Total 



276.1 ±23.5 



the scattered neutron but the cross sections are not known 
precisely. A ^^°Po^^C source was employed to study the 
^^C(Q;,n)^^0 reaction and tune a simulation using the cross 
sections from Ref. [10, 11]. We find that the cross sections for 
the excited 1*^0 states from Ref. [10] agree with the ^^Opo^^C 
data after scaling the P' excited state by 0.6; the 2"'' excited 
state requires no scahng. For the ground-state we use the cross 
section from Ref. [11] and scale by 1.05. Including the ^^''Po 
decay-rate, we assign an uncertainty of 1 1 % for the ground- 
state and 20% for the excited states. Accounting for e(-Ep), 
there should be 182.0 ±21.7 ^'^Cia^nf^O events in the data. 

To mitigate background arising from the cosmogenic beta 
delayed-neutron emitters ^Li and '^He, we apply a 2 s veto 
within a 3-m-radius cylinder around well-identified muon 
tracks passing through the LS. For muons that either deposit 
a large amount of energy or cannot be tracked, we apply a 2 s 
veto of the full detector. We estimate that 13.6 ± 1.0 events 
from ^Li/^He decays remain by fitting the time distribution of 
identified ^Li/^He since the prior muons. Spallation-produced 
neutrons are suppressed with a 2 ms full-volume veto after a 
detected muon. Some neutrons are produced by muons that 
are undetected by the OD or miss the OD but interact in the 
nearby rock. These neutrons can scatter and capture in the LS, 
mimicking the Fe signal. We also expect background events 
from atmospheric neutrinos. The energy spectrum of these 
backgrounds is assumed to be flat to at least 30 MeV based on 
a simulation following [12]. The atmospheric v spectrum [13] 
and interactions were modeled using NUANCE [14]. We ex- 
pect fewer than 9 neutron and atmospheric v events in the 
data-set. We observe 15 events in the energy range 8.5 - 
30 MeV, consistent with the Umit reported previously [15]. 

The accidental coincidence background above 0.9 MeV is 
measured with a lO-ms-to-20-s delayed-coincidence window 
to be 80.5 ±0.1 events. Other backgrounds from (7,n) inter- 
actions and spontaneous fission are negligible. 

Anti-neutrinos produced in the decay chains of ^■^^Th and 
^■^^U in the Earth's interior are limited to prompt ener- 
gies below 2.6 MeV. The expected geo-neutrino flux at the 
KamLAND location is estimated with a geological reference 
model [9], which assumes a radiogenic heat production rate 
of 16 TW from the U and Th-decay chains. The calculated 
fluxes for U and Th-decay, including a suppression factor of 
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FIG. 1: Prompt event energy spectrum of 77e candidate events. 
All histograms corresponding to reactor spectra and expected back- 
grounds incorporate the energy-dependent selection efficiency (top 
panel). The shaded background and geo-neutrino histograms are cu- 
mulative. Statistical uncertainties are shown for the data; the band on 
the blue histogram indicates the event rate systematic uncertainty. 



0.57 due to neutrino oscillation, are 2.24x 10^ cm ^s ^ (56.6 
events) and 1.90x lO^cm^-^s^^ (13.1 events), respectively. 

With no Fe disappearance, we expect 2179 ± 89 (syst) 
events from reactors. The backgrounds in the reactor energy 
region listed in Table II sum to 276.1 ±23.5; we also expect 
geo-neutrinos. We observe 1609 events. 

Figure 1 shows the prompt energy spectrum of selected 
Vf, events and the fitted backgrounds. The unbinned data is 
assessed with a maximum likelihood fit to two-flavor neu- 
trino oscillation (with 6*13 = 0), simultaneously fitting the geo- 
neutrino contribution. The method incorporates the abso- 
lute time of the event and accounts for time variations in 
the reactor flux. Earth-matter oscillation effects are included. 
The best-fit is shown in Fig. 1. The joint confidence in- 
tervals give Am2^ = 7.58l°j^(stat)to;i5(syst) x IQ-^eV^ 
and tan^ 6I12 = 0.56l|J;J?(stat)+°:J°(syst) for tan^ 6ii2<l. A 
scaled reactor spectrum with no distortion from neutrino os- 
cillation is excluded at more than 5(7. An independent anal- 



5+0.22 
-0.20 



ysis using cuts similar to Ref. [2] gives Amji = 7.66 
10-5 eV2 and tan^ d^^ = 0.52ta \t 

The allowed contours in the neutrino oscillation parame- 
ter space, including Ax^-profiles, are shown in Fig. 2. Only 
the so-called LMA-I region remains, while other regions pre- 
viously allowed by KamLAND at ^2.2a are disfavored at 
more than 4(t. For three-neutrino oscillation, the data give 
the same result for Am^i, but a slightly larger uncertainty on 
9i2. Incorporating the results of SNO [16] and solar flux ex- 
periments [17] in a two-neutrino analysis with KamLAND as- 
suming CPT invariance, gives Am'^i = T.SOlg 21 ^ lO^^ eV^ 
and tan^ 6I12 = 0.47t|5;!^^. 

To determine the number of geo-neutrinos, we fit the nor- 
malization of the Ve energy spectrum from the U and Th- 
decay chains simultaneously with the neutrino oscillation pa- 
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FIG. 2: Allowed region for neutrino oscillation parameters from 
KamLAND and solar neutrino experiments. The side-panels show 
the Ax'^-profiles for KamLAND (dashed) and solar experiments 
(dotted) individually, as well as the combination of the two (solid). 
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FIG. 3: Ratio of the background and geo-neutrino-subtracted 77e 
spectrum to the expectation for no-oscillation as a function of 
Lq/E. I/O is the effective baseline taken as a flux-weighted aver- 
age (Lo = 180 km). The energy bins are equal probability bins of the 
best-fit including all backgrounds (see Fig. 1). The histogram and 
curve show the expectation accounting for the distances to the indi- 
vidual reactors, time-dependent tlux variations and efficiencies. The 
error bars are statistical only and do not include, for example, corre- 
lated systematic uncertainties in the energy scale. 



rameters using the KamLAND and solar data. There is a 
strong anti-correlation between the U and Th-decay chain 
geo-neutrinos and an unconstrained fit of the individual con- 
tributions does not give meaningful results. Fixing the Th/U 
mass ratio to 3.9 from planetary data [18], we obtain a 
combined U-nTh best-fit value of (4.4± 1.6)x lO'^cm^^s^^ 
(73 ± 27 events), in agreement with the reference model. 

The KamLAND data, together with the solar v data, set an 
upper limit of 6.2 TW (90% C.L.) for a reactor source at 
the Earth's center [19], assuming that the reactor produces a 
spectrum identical to that of a slow neutron artificial reactor. 

The ratio of the background-subtracted candidate events, 
including the subtraction of geo-neutrinos, to no-oscillation 
expectation is plotted in Fig. 3 as a function of Lq/E. The 
spectrum indicates almost two cycles of the periodic feature 
expected from neutrino oscillation. 

In conclusion, KamLAND confirms neutrino oscillation, 
providing the most precise value of Am^i to date and im- 
proving the precision of tan^ 9i2 in combination with solar v 
data. The indication of an excess of low-energy anti-neutrinos 
consistent with an interpretation as geo-neutrinos persists. 
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